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We study the interaction between electron and acoustic phonon in a Rashba spin-orbit coupled 
two dimensional electron gas using Boltzmann transport theory. Both deformation potential and 
piezoelectric scattering mechanisms are considered in the Bloch-Gruneisen (BG) as well as in the 
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dence of resistivity in the BG and EP regimes has been discussed. We find effective exponent of the 
temperature dependence of the resistivity in the BG regime decreases due to spin-orbit coupling. 
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I. INTRODUCTION 

With the rise of the promising field of spintronics, two 
dimensional electron gas (2DEG) with spin-orbit interac- 
tion (SOI) in semiconductor heterostructures has drawn 
much attention in recent years 1 - - — . The importance of 
this field was first realized when Datta and Das gave a 
proposal of Spin Field Effect Transistor 4 . Several stud- 
ies on spintronics have been performed in recent years 
from both theoretical and experimental viewpoints. One 
major aim is to manipulate the spin degree of freedom 
of charge carriers in semiconductor nanostructures^ so 
that spin-based device technology^ and quantum infor- 
mation processing^ can be developed in near future. The 
SOI is an intrinsic phenomena present in semiconduc- 
tors. Mainly, there are two kinds of SOI present in the 
semiconductor heterostructures we come across in the lit- 
erature. One of them is the Rashba spin-orbit interaction 
(RSOI)^ which originates from the inversion asymmetry 
of the confining potential in semiconductor heterostruc- 
tures. The RSOI is proportional to the magnitude of the 
electric field internally generated due to the band bend- 
ing or externally applied gate voltage. It can be tuned 
by applying a gate voltag o 9 i . Another kind of SOI is 
the Dresselhaus spin-orbit interaction^ which originates 
from the bulk inversion asymmetry of the host crystal. It 
entirely depends on crystal property and it is not tunable. 

Various electronic and transport properties of a 2DEG 
will be modified in the presence of the SOI. We mainly fo- 
cus on the RSOI in 2DEG systems such as AlGaAs/GaAs 
heterostructure. The coupling between electron and 
phonon becomes stronger when the SO coupling constant 
is large enough and due to this fact the effective mass is 
also increased^. The critical temperature of supercon- 
ductors can be controlled by RSOI when the Fermi energy 
is small as compared to the characteristic energy scale of 
RSOI 13 . The RSOI can produce infinite number of bound 
states in a 2DEG with short-range impurity potentials 1 ^. 
At low frequency in the presence of RSOI the universal- 
ity of spin Hall conductivity can be broken by consid- 
ering the contribution of electron-phonon interaction to 
the spin-vertex correction 1 ^. The mobility, polaron mass 



correction and polaron binding energy can be changed 
significantly due to the presence of the RSOI 1 -' 17 . The 
relaxation time for various impurity potentials of low- 
dimensional semiconductor structures with SOI has been 
studied^. 

The interaction between electron and phonon plays a 
very crucial role in determining the transport properties 
of a 2DEG and it has a finite contribution to the mo- 
mentum relaxation time of the charge carriers. Other 
contributions 1 ^ - — come from disorders, impurities, etc. 
There are three distinct temperature regimes: a) Bloch- 
Gruneisen (BG), b) Equipartition (EP) and c) Inelas- 
tic phonon scattering. The BG temperature Tbg can 
be defined^! by the relation ksTsc = 2hkpv s , where 
ks is the Boltzmann constant, v s is the phonon veloc- 
ity and kp is the Fermi wave vector. For typical elec- 
tron density in e ~ 10 15 m~ 2 ) in a 2DEG system, Tbg 
is around 6.2 K. In the BG regime, a direct manifesta- 
tion of the acoustic phonon-dominated transport prop- 
erty is the strong change in the temperature depen- 
dence of the resistivity. The existence of the BG regime 
has been confirmed for 2DEG experimentally 24 . The 
problem of electron-phonon interaction in a 2DEG con- 
fined in semiconductor heterostructures has been studied 
extensively 2 ^—. Recently, the phonon-dominated trans- 
port properties in the BG regime have been studied in 
graphene both theoretically 3 ^ - — and experimentally 3 ^. 

In the present work we would like to investigate the 
influence of the RSOI on momentum relaxation time due 
to electron-phonon interaction and hence on the trans- 
port properties of 2DEG systems. We consider both the 
cases for perfect 2DEG and quasi-2DEG (usually found 
in semiconductor heterostructures). In the former case 
2D phonon wave vector q couples with 2D electron wave 
vector k and in the latter case the coupling between 3D 
bulk phonon wave vector Q = (q, q z ) and 2D electron 
wave vector k is considered. We consider two mechanisms 
of electron-phonon interaction, namely deformation po- 
tential (DP) and piezoelectric (PE) potential scattering 
separately. In all the cases linear temperature depen- 
dence of inverse relaxation time (IRT) is found in the 
high temperature (EP) regime. But in the BG regime 
for perfect 2DEG and quasi-2DEG we find analytically 
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resistivity is proportional to T 4 and T 5 , respectively, in 
the case of DP scattering mechanism. On the other hand, 
p ~ T 3 for PE scattering mechanism in a quasi-2DEG. 
Our numerical calculations reveal that this exponent of 
T strongly depends on density and the SO coupling con- 
stant. In fact, the exponent of T decreases due to pres- 
ence of spin-orbit coupling. We also discuss the resistivity 
as a function of the SO coupling constant in the BG and 
EP regimes. 

This paper is organized as follows. In section II we 
derive all the theoretical results and discuss all the nu- 
merical analysis for perfect 2DEG. In section III both 
DP and PE scattering mechanisms have been taken into 
account for quasi-2DEG and we discuss analytical and 
numerical results in detail. We summarize all the results 
in section IV. 



H cp — _DV • u(r), where D is the deformation-potential 
coupling constant and the lattice displacement vector 
u(r) is given by 
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Here, w q = v s q is the phonon frequency with the wave 
vector q and the sound velocity v s , a q and a q are phonon 
creation and annihilation operators, respectively. Also, 
e q is a unit vector in the direction of the phonon polar- 
ization. The corresponding Hamiltonian can be written 
as 



H cp (r) = £ [C. 
q 



q"q>- 
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II. ELECTRON-PHONON SCATTERING IN A 
PERFECT 2DEG 

A. Theoretical model 

We consider a 2DEG with the RSOI in the xy plane 
in a semiconductor heterostructure. The single-particle 
Hamiltonian of this system is given by 



P a i 



(1) 



where p is the two-dimensional momentum operator, m* 
is the effective mass of an electron, ao is the unit 2x2 
matrix, a is the Rashba spin-orbit coupling constant and 
o~ x ( y ) are the Pauli matrices. The eigenenergies are given 
by 



£ A(fc) = ^7 + Aa|k|, 
with the corresponding normalized eigenspinors 



(2) 



(3) 



Here, A = ± represents the upper and lower en- 
ergy branches, A is the area of the system and <j> = 
tan -1 (k x I 'k y ) . The density of states for the two energy 
branches are given by-i£ 



D ± (e) = ^(1 I 



V e + £a > 

+ D J-^e(-e)Q(e + e a ). (4) 

Here, Dq = m*/(nh 2 ), G(x) is the unit step function 
and e a = m*a 2 /(2h 2 ) is the characteristic energy scale 
of RSOI. 

The Hamiltonian for electron-phonon interaction in the 
case of deformation potential coupling can be written as 



where C q = D v / h/2MNujq(ie q ■ q). 

The energy dependent relaxation time for electrons in 
a given energy branch A can be written as 



1 



r x (e) 
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where 9^ is the scattering angle between the two mo- 
mentum vectors k and k', e\ is given by Eq. |2]), P££, 
is the transition rate for scattering of an electron from 
a state |k,A) to |k',A') and /(e) = [e^"^ + l]" 1 is 
the Fermi-Dirac distribution function with /3 — l/(fcsT). 
The chemical potential [i at finite temperature T can be 
obtained self-consistently from the following normaliza- 
tion condition^: 



/■oo />0 

J deD a f{e) + J deDv 



e + e c 



(8) 



where n e is the electron density. At T — 0, the above 
equation reduces to ep — ep — 2e Q , where tp and 
is the Fermi energy of a 2DEG in presence and absence 
of RSOI, respectively. Thus, the reduction in the Fermi 
energy due to RSOI is 2e Q . 

We consider the interaction between electron and 
acoustic phonon due to deformation potential coupling. 
The transition rate due to electron-phonon interaction 
can be written as 



pAA' 
Mck' 



2tt 

T 



E 

q 



C. 



AA' \2 



(A^ q + l)5(e' x , -e x + fiw q )l , 



(9) 



where is the matrix element for the acoustic phonon 
and is given by 

I^aa',2 £> 2 % l + AA'cosfl 

|Cq 1 = 2A^T 2 5xX '- (10) 

Here, 6 = 0^*, D is the deformation potential coupling 
constant and p a = NM/A is the mass per unit area. The 
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appearance of 5\\i is due to the fact that the electron- 
phonon interaction given by Eq. is spin- independent. 
Also, Nq = [exp(f3hajq) — is the phonon occupation 
number. The first and second terms on the right hand 
side of Eq. ^ correspond to the absorption and emission 
of a phonon with energy hujq, respectively. Within the 
small-angle scattering approximation (q = 2k sin(0/2)), 
the matrix element for intra-branch scattering (A = A') 

becomes |C^ V | 2 = 2 ^ p hq v (1 - f^). The similar matrix 
element is obtained for a single layer graphene 3 - 3 -. 

Before presenting the numerical results, we present 
how IRT depends on T in the EP and the BG regimes. 
At high temperature (EP regime), the phonon energy 
is much smaller than the thermal energy i.e. hujq <C 
fcflT. We neglect hujq term in the delta functions i.e. 
8(e' — e =F hujq) ~ 5(e' — e). Again, in this temperature 
limit the Bose occupation factor A^ q can be approximated 



as N a 



1 ~ kBT/hujq. When ep > e a , the total 



relaxation time at high temperature is 



1 



m*D 2 



k B T. 



(11) 



Therefore, in the EP regime, the IRT depends linearly 
on temperature. 

Now we want to see how resistivity depends on tem- 
perature and a at low temperature (BG regime) where 
hujq ~ ksT . In the BG regime, the IRT strongly de- 
creased because the phonon population decreases expo- 
nentially for phonon absorption and the sharp Fermi dis- 
tribution prohibits phonon emission. To see the temper- 
ature dependence of the resistivity in the BG regime, it 
is convenient to calculate IRT averaged over energy, as 
used for graphene^ 3 ., which is given by 



1 \ _ Jjeg^hgg 
f d eD x (e)[-SM] 



(12) 



Therefore, the resistivity of a 2DEG with RSOI can be 
calculated from the following equation: 



n e e^ 



(13) 



where (1/r) = £ A (l/r A ). _ 

When the temperature is very low we can make the 
following approximations: (i) the phonon energy is com- 
parable with the thermal energy i.e ksT < fvjjq < ej? 
and /(e) [1 - /(e ± hwq)} ~ fiw q (iV q + 1/2 ± 1/2) J(e - ep). 
After taking these approximations, we obtain 



2AD 
k H T 



it 



e F + e c 



X / dO(l-COS9)\Cq\ 2 LOqNq(Nq + l). 



(14) 

We can now convert the integration over 9 into q by us- 
ing the relation (based on small-angle scattering) q = 



2k F sin | . Substituting this and after a straightforward 
calculation we finally obtain 



D, 
2 
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where k F — \j fc 2 + k° F 2 ^k a , k a = m*a/h 2 is the Rashba 
wave vector, k F = ^/2ixn e is the Fermi wave vector of a 
2DEG without RSOI and C(4) = vr 4 /(90). In this tem- 
perature regime energy averaged IRT is proportional to 
T 4 and hence p ~ T 4 . The T 4 scaling law is also found 
in other perfect 2D system such as graphene^. As we 
will see in the numerical calculations, the exponent of 
the temperature dependence of the resistivity is reduced 
due to SOI coupling. 



B. Numerical results 

In this section we determine the IRT, resistivity and 
their dependence on energy, temperature, SOI coupling 
constant etc. We solve Eq. ([7} numerically. For our 
numerical calculation, we use m* = 0.067mo with mo 
is the free electron mass, a = ao = 1CP 11 eV-m, v s — 
5.3 x 10 3 ms _1 and the electron density n = 10 15 m~ 2 . 

In Fig. 1 we have plotted IRT as a function of energy 
at a fixed temperature T = 1 K for densities n e = no and 
n e = 5uq. It could be seen from Fig. 1 that there is a dip 
in the IRT. The dip occurs due to the sharpness of the 
Fermi distribution function at low enough temperature. 
It is clear from Fig. 1 that the dip occurs exactly at 
e = €° F when a = 0. As a increases the position of the 
dip appearing at energies lower than e° F . The dip occurs 
when e — p, where p is the chemical potential which we 
evaluate numerically by solving Eq. ([5]) . For the electron 
density no, the Fermi energy without RSOI is e^=3.6 
meV. When a = 0, u = e F = 3.6 mcV. When a = 3ao, 
(f = 2.8 meV and p — 0.8e^ so the dip occurs exactly 
at e = 0.8e^. Similarly, for a = 5ao, ep = 1.4 meV and 
u = 0Ae° F and consequently the dip occurs exactly at 
e = 0Ae° F . The IRT of the 2DEG system with RSOI is 
reduced compared to the absence of SOI. 

The effective exponent (y) of the temperature depen- 
dence of the resistivity strongly depends on the density 
and the RSOI coupling constant. We estimate the ef- 
fective exponent from the log-log plot of the resistivity 
versus T in Fig. 2 for a — and a = ao with differ- 
ent densities. In the BG regime (T ~ 1 — 3 K), we find 
v = 3.763,4.111 and 4.235 for n e = 3no,5no and 7no, 
respectively, when a = 0. On the other hand, we find 
v = 2.829,3.203 and 3.433 for n e = 3no,5no and 7no, 
respectively, when a = ao- The exponent v is increasing 
with electron density. 

In the high temperature limit around T — (18 — 40) 
K, we get v = 0.824,0.998 and 1.019 for densities n e = 
3no,5no and 7no, respectively. The v is also increasing 
with n e slowly and the equipartition result p ~ T is re- 



FIG. 1: (Colo r online) Plots of the IRT [in units of l/r = 
2m* D 2 \/2~imo / (nh 2 p a v s )] versus energy for different values 
of a. The solid, dashed and dot-dashed lines correspond to 
the a — 0, a — 3ao and a = 5a?o, respectively. For better 
visualization, the solid lines in the left and the right panels 
have been reduced by a factor of 3 and 4, respectively. 



T (K) T (K) 



FIG. 3: (Color online) Plots of the resistivity with tempera- 
ture at a fixed density for different values of a. Here, solid 
and dashed lines correspond to the a = and a = ao, respec- 
tively. 




covered. The exponent v does not change due to a in the 
high temperature limit. 

In Fig. 3, we plot p versus T at a fixed density for 
different values of a. In the each panel of Fig. 3, we 
consider a — and a — ao cases for densities n e = 3no 
and n e = 5no- Figure 3 clearly shows that the slope of the 
curve for a = is greater than that for a = olq case in the 
BG regime. When n e — 3rto, we estimate v — 3.763 and 
v = 2.829 for a — and a = ao, respectively. Similarly, 
when n e — 5n v — 4.111 and v = 3.203 for a — 
and a = ao, respectively, at very low temperature. The 
values of v differ significantly between a = and finite 



FIG. 4: (Color online) Plots of p of a perfect 2DEG versus a 
for two fixed densities. Here solid and dashed lines correspond 
to n e = 3no and n e = 7no, respectively. 



a at very low temperature. 

In Fig. 4, we plot resistivity of a perfect 2DEG versus 
a for two different densities. In both the regimes, p in- 
creases with a. But the rate of increase of p is high for 
low electron density compared to that of the high density. 
In the BG regime, the resistivity is saturated at high a. 
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III. ELECTRON-PHONON SCATTERING IN A 
QUASI-2DEG 

In this section we consider three dimensional bulk 
phonon with wave vector Q = (q, q z ) interacts with the 
two dimensional electron wave vector k = (k x ,k y ). The 
component of Q in the xy plane obeys the conserva- 
tion of momentum i.e k — k' = q and the z-component 
must be integrated out. In semiconductor heterostruc- 
ture the electrons move in the x-y plane in presence of 
a triangular potential in the z direction. It is also as- 
sumed that the lowest energy level is occupied by the 
electrons. Generally the wave function can be written as 
tjj(r) = ip(x,y)C,o(z). Here the variational wave function 
Co(z) is given by, Co( z ) = \/b 3 j2ze~ hz l 2 . The variational 
parameter— b is given by b = (487rm*e 2 / 'EoKoti 2 ) 1 / 3 (n c + 

x 1/3 

lln e /32 J , where kq — 12.9 is the static dielectric con- 
stant of GaAs, £o is the free space permittivity and n c 
is the depletion charge density in the channel. In this 
section we will also discuss the piezoelectric scattering 
along with the deformation potential scattering. In this 
case the IRT can be written as 

x J dq z \I(q z )\ 2 \C^'f{N Q 5(e' x , - e A - hw Q ) 
+ (N Q + l)5(e' y -e x + hw Q )} ' (16) 



where the term |/(g^)| 2 , called the form factor, is re- 
sponsible for getting transition rate from three dimen- 
sional bulk phonon state to two dimension. The exact 
form of this term for a triangular potential is given by 



\I(<h 



Jdz( 2 (z)e l i* z \ 2 =b 6 /(b 2 + q 



2\3 



DP scattering as given by 

1 3 m*bD 2 



Tie 



32 p m h 3 v\ 



k B T. 



(19) 



On the other hand, in the BG regime, we obtain the 
following expression for the energy averaged IRT (see ap- 
pendix Al): 



1 



Do 



DP 



it 



D 2 5\((5)(k B T) 5 



- e a J p m v s (v s h) 5 (kp) 



where C(5) = 1.037. 



(20) 
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FIG. 5: (Color online) Plots of the resistivity of a quasi-2DEG 
due to DP scattering versus temperature for different values 
of the density. Here, solid, dotted and dashed lines represent 
n e = 3no, n e = 5no and n e = 7no, respectively. 



A. Deformation potential scattering 

The matrix element in this case is given by 



where p rn is the mass density. At high temperature, we 
have Nq ~ Nq + 1 ~ ksT '/ (hv s Q). Inserting this matrix 
element into Eq. (fir!)) , we get 



2m * D 'l T J-^) kB T 



r ± (e) n 2 h 3 p m v] 

x / draVl-i 2 / dq z \I(q z )\ 2 , (18) 
Jo Jo 

with x — q/2k. Performing the integrations over x and 
q z , we finally obtain the total IRT in the EP regime for 



B. Piezoelectric scattering 

Piezoelectricity is nothing but generation of polariza- 
tion due to the application of a strain to a crystal without 
inversion symmetry. Due to lattice vibration a potential 
can be generated in such crystals and electrons are scat- 
tered by this kind of potential. To calculate the IRT due 
to the PE scattering, we can use Eq. (fT6|). Following 
Refi2£, the matrix elements of the Rashba system are 
obtained as 



pe,\,\',2_ {e-hii) 2 h 1 + AA' cos6> 



A l{t) (q,q z ), (21) 



where Ai(q,q z ) = 9q 2 z q 4 /2(q 2 z + q 2 ) 3 and A t {q,q z ) = 
(8g^ 2 + g 6 )/4(g 2 + g 2 ) 3 . In Eq. (J2TJ) the value PE tensor 



component hn is 1.2 x 10 V/m and v s ;( t ) is the lon- 
gitudinal (transverse) component of sound velocity. By 
inserting this matrix element into Eq. (|16p and doing 
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and 
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FIG. 6: (Color online) Plots of the resistivity of a quasi-2DEG 
due to DP scattering versus a for various densities. Here, 
solid and dashed lines correspond to n e = 3no and n e = 7no, 
respectively. 



some straightforward calculations in the high tempera- 
ture regime we get, 



9 to* (e/iu) 2 



and 



T ±(e) 8tt 2 £3 p m v 2 sl 

d6 sin 2 9Fi{q) 

1 to* {ehi^Y 



it 



(22) 



it 



r*(e) ' " 16tt 2 ft 3 p m v 2 t 

x f d9 sin 2 d F t (q). 
Jo 

Here i 7 / (q) and Ft (?) are respectively given by 



k B T 



(23) 
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2 g;g 
(<z 2 + <7 



7r 1 + 6k+12k 



2 ' 2k 3 



and 



16q 



dg*|J(g z )| 



(1 + k) c 



2 8g 4 <7 2 + g 6 



(24) 



(^Tq 2 ) 4 16« 
13 + 78k T 72k 2 + 82k 3 + 36k 4 + 6k 5 
(1 + k) 6 



,(25) 



where k = q/b. Now, we are assuming that the quasi- 
2DEG is very thin, i.e k« 1. So Fi(q) and Ft{q) can be 
approximated as Fi(q) ~ n/16q and i^t(g) — ISir/Wq. 

Substituting Fi(q), F t (q) and q — 2ksin(9/2) in Eqs. 
([22| and (|23|) and integrating over 0, we obtain 



3 to* (e/m) 2 1 
32 7rft 3 p m w 2 , /c± 



fcsT (26) 



1 13 to* (e/ii 4 ) 2 1 

7f(e) ~ 192 irh 3 Pm v 2 st fc± 



fc B T. (27) 



The total IRT for longitudinal and transverse cases are 
given as 



1 



and 



3 to* (eh 14 ) 2 k B T 

Ti(e) ~ IQirh? p m v 2 sl k° F 



1 13 to* (ehu) 2 k B T 




T f (e) 96 7rft 3 p m v 2 t k% V e + e Q 



(28) 



(29) 



The total IRT in the PE scattering case can be written as 
1/t pe = l/rj+2/rt. In the high temperature regime, the 
IRT is proportional to k B T and inversely proportional to 

In the low temperature regime we calculate IRT av- 
eraged over energy. The detail calculations are given in 
appendix A2. In this case, we have following expressions 



45 to* (eh 14 ) 2 3!C(3) (k B T\?> 



PE 
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kp 



(30) 



and 



PE 



59 to* (eh 14 ) 2 3!C(3) (k B T\t 



1024 irh p m (hv st y \ k$ J 



Here, C(3) = 1.202. Equations (HP and ([21]) show that 
the energy averaged IRT due to PE scattering is propor- 
tional to T 3 and hence the resistivity is also proportional 
to T 3 . 



C. Numerical results 

In this section we discuss resistivity due to both DP 
and PE scattering mechanisms. We solve Eq. (fT7J|) nu- 
merically using the matrix elements for DP and PE scat- 
tering given in Eqs. (fl~7l) and (|2"Tj) . respectively. For 
the numerical calculation we set v s i = 5.31 x 10 3 ms _1 , 
v st = 3.04xl0 3 ms-\ n c = 5xl0 14 m~ 2 , p m = 5.12xl0 3 
Kgm -3 and D = 12 eV. The other parameters are the 
same as given in Sec. II. 

In Fig. 5, we show the temperature dependence of the 
resistivity of the quasi-2DEG due to the DP scattering for 
fixed values of a — and a — a$ with different densities. 
In the BG regime with a = 0, the effective exponents 
of T are v = 4.459,4.887 and 5.085 for densities n e = 



7 




4 5 6 78910 
T(K) 




T(K) 



FIG. 7: (Color online) Plots of the resistivity due to PE 
scattering versus temperature for different values of the den- 
sity. Here, solid, dotted and dashed lines represent n e = 3no, 
n e = 5no and n e = 7no, respectively. 



BG regime. 

In Fig. 7, we plot the temperature dependence of the 
resistivity of the quasi-2DEG due to the PE scattering for 
fixed values of a — and a — ao with different densities. 
When a = 0, the exponents are v — 2.537, 2.830 and 
3.002 for densities n e — 3no,5no and 7uq, respectively, 
in the BG regime. Similarly, the effective exponents are 
v = 1.962,2.208 and 2.398 for densities n e = 3no,5no 
and 7no, respectively, when a — ao- It clearly shows 
that the exponent v is decreased due to the presence of 
the RSOI. 

For PE scattering case we also calculate the resistivity 
of a quasi-2DEG as a function of a by solving Eq. (fTB|) 
and plot the results in Fig. 8. Figure 8 shows that p 
increases with a in both the regimes. The other features 
are similar to the DP scattering mechanism. 

Comparing Fig. 6 and Fig. 8, one could see that the 
resistivity due to DP and PE potential are in the same 
order in the BG regime but p due to DP is dominating 
over PE in the EP regime. 



IV. SUMMARY 
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FIG. 8: (Color online) Plots of the resistivity of a quasi-2DEG 
due to piezoelectric scattering versus a for various densities. 
Here, solid and dashed lines correspond to n e = 3no and 
n e = 7no, respectively. 



In this work we have investigated the effect of the RSOI 
on the momentum relaxation time due to the electron- 
phonon scattering in a 2DEG. We have considered both 
perfect 2DEG and quasi-2DEG. We have also considered 
both the DP and PE scattering mechanisms responsi- 
ble for the electron-phonon interaction separately. The 
temperature dependence of the resistivity has been cal- 
culated in both EP and BG regimes. We have found 
through approximate calculations that the resistivity of 
a perfect 2DEG is proportional to T 4 in the BG regime 
for the DP scattering. On the other hand, p ~ T 5 for 
DP scattering and p ~ T 3 for PE scattering in a quasi- 
2DEG. We have also recovered the linear temperature 
dependence of the resistivity in the EP regime of all the 
cases. Our numerical analysis showed that the effective 
exponent {v) of the temperature dependence of the resis- 
tivity strongly depends on the electron density n e and SO 
coupling constant a. For the DP and PE scattering, the 
values of the effective exponent of T in the BG regime 
for different values of n e and a are summarized in Table 
I. 



3no , 5no and 7no , respectively. Similarly, for a = ao ■ 
we get the exponents as v = 3.6499,4.058 and 4.359 for 
densities n e = 3no,5no and 7no, respectively. In this 
case, v is greater than the perfect 2DEG case discussed 
in the previous section. This increase in v is due to the 
finite thickness of quasi-2DEG in the z direction. In the 
EP regime, the v is closed to one as we increase n e . 

In Fig. 6, we plot the resistivity of a quasi-2DEG due 
to DP scattering as a function of a at fixed temperature 
and density. Figure 6 shows that the p increases rapidly 
with a in both the regimes but it is faster for low electron 
density and gets saturated after certain value of a in the 





v for 2DEG 


v for quasi-2DEG 


Density 


DP 


DP 


PE 


(n e ) 


a = 


a = oto 


a = 


a = Qo 


a = 


a = ao 


3no 


3.763 


2.829 


4.459 


3.649 


2.537 


1.962 


5no 


4.111 


3.203 


4.887 


4.058 


2.830 


2.208 


7n 


4.235 


3.433 


5.085 


4.359 


3.002 


2.398 



TABLE I: The effective exponent of the temperature depen- 
dence of the resistivity in the BG regime for various values of 
n e and a. 

The variation of p with a has also been discussed for 
all the cases. It is found that the p increases with a in 



s 



both the regimes. The rate of increase is faster in low Using the result J de p e™ Aq^q+I) = n\((n)(kBT) n+1 
electron density case. with is the Riemann zeta function, we finally obtain 



Appendix A 

In this appendix, we derive energy-averaged IRT due to 
DP and PE scattering in a quasi-2DEG in the BG regime. 
At very low temperature, following the approximations 
used in section 11(A) and using Eq. (fT6|) . we get 



1 



m* / 
^3 V 



it 



i 



e a J k B T 



dO(l - cos( 



x J dq z \I{q z )\ 2 \C^ q fhu: Q N Q {N Q + l). (Al) 

When T is very low, the phonon wave vector q is very very 
small compared to the Fermi wave vector i.e. q <C 2k F - 
For very thin quasi-2DEG, |/(g z )| 2 can be approximated 
as |J(g z )| 2 ~ 1 since q z <C b. With all the approxima- 
tions taken into account, Eq. (|A1|) can be approximated 
further as 



1 



2tt 2 ^ 3 k 



it 



i 



knT 



J dqdq z q 2 \C^ q fhLO Q N Q {N Q + 1). (A2) 



1. Deformation Potential scattering 

Inserting the matrix element given in Eq. (|17|) into Eq. 
A2]) . we get 



1 



lV 2 k$ 3 i 1 T \le° F 6 -eJ k B T 



X [ dqdq z q 2 {hw Q ) 2 N Q {N Q + l). (A3) 



T= / DP 



Since the phonon dispersion relation is e p = 
hv s q 2 + q 2 , we can make the following transformation: 
dqdq z — > e p de p d(f>/(hv s ) 2 with q — e p cos (f)/(Hv s ) and 
Qz = e P sm(j)/(hv s ). With these transformations, Eq. 
(IA2I) reduces to 



1 



m*D 2 



r± / dp 4ir 2 h 3 p m v 2 kp 3 

X ~kbr I de p4 N Q( N Q + 1 



(A4) 



1 \ D °,> 
^)dp^- [1T 



D 2 5!C(5)(fc B T) 5 



p m v s (v s h) 5 (kp) 3 



(A5) 



2. Piezoelectric scattering 



Using the matrix elements for PE scattering given in 
Eq. (f2"Tj) and Eq. (|A2|) . we obtain the energy-averaged 
IRT for longitudinal and transverse cases, respectively, 



1 



9 m*(ehu) 



tt I PE 



it 



8ir 2 p m h(Hvsi) 4 n, F 
J de p e p N Q (N Q + 1) J cos 6 </>sin 2 



1 



knT 



and 



1 m*(ehi4) 2 



it 



8tt 2 Pm h{hv st Ykf 
x / (cos 8 (j> + 8 cos 4 (j> sin 4 



1 



x de p e p N Q (N Q + l). 



£ F - tot' ksT 



(A6) 



After doing the integration over e p and <p, we finally ob- 
tain the following expressions 



45 m* (ehu) 2 3!<C(3) /k B T\ 3 
512 ^ft p m (hv sl ) 4 \~kYJ 



if I PL 



and 



59 m* (e/ii 4 ) 2 3!C(3) /fc B T\3 
rf/pE ~ 1024^ p m (frv st ) 4 \~k F r 



(A8) 
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